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Abstract 


Surface incompressibility, also called iuextensibility, imposes a zero-surface-divergence con¬ 
straint on the velocity of a closed deformable material surface. The well-posedness of the 
mechanical problem under such constraint depends on an inf-sup or stability condition for 
which an elementary proof is provided. The result is also shown to hold in combination with 
the additional constraint of preserving the enclosed volume, or isochoricity. These continu¬ 
ous results are then applied to prove a modified discrete inf-sup condition that is crucial for 
the convergence of stabilized finite element methods. 

Keywords: Inextensibility, Surface incompressibility, Volume preservation, Inf-sup 
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i 1. Introduction 

The rate of change of the volume |cu| of an infinitesimal piece of continuum u located at 
point x at time t is given by [l| 



( 1 ) 


There exist materials that preserve volume exactly, or within experimental accuracy, and are 
modeled as incompressible. The incompressibility constraint (V • u = 0) materializes in the 
equations of motion as a reaction force, which is the gradient of an unknown pressure field 
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p. The pressure field is a uniquely determined element of Q = L 2 (Q) because the diver gen ce 


operator V- is surjective onto Q. The surjectivity is equivalent to the inf-sup condition 

v 


inf 


sup . 

qsQ veH^O) IMIl ||9||o 

If the continuum is restricted at the boundary, more precisely, if the velocity field normal to 
the boundary is constrained all over dfl, then the mean value of the pressure is undetermined. 
This reflects mathematically in that to satisfy (J2]) the pressure space Q must be chosen as 
a strict subspace of L 2 (fl), for example Lq(Q) which consists of functions with zero mean. 

If we now consider a smooth, closed, orientable surface T evolving in R 3 , the analogous 
property to volume preservation is area preservation, also called inextensibility or surface 
incompressibility, which is indeed exhibited by many relevant materials. Our main inter¬ 
est is in lipid membranes [i, 5, 6, 7], which are area-preserving surface fluids |8f, but the 
inextensibility constraint can also hold in other, fluid or solid, material surfaces. 

The rate of change of the area \w\ of an infinitesimal piece of surface zu located at x at 
time t and moving along a velocity field u is given by 

d\w 


> 0 . 


, 3 ] 

( 2 ) 


dt 


= \w\ Vr • u(x, t ) 


( 3 ) 


where we have introduced the surface divergence operator Vr - - The velocity field of an 
inextensible surface T must thus satisfy 


Vr • u(x, t) = 0 a.e. in T . 


( 4 ) 


Remark 1. We adopt here the operators of tangential calculus as presented by D elf our and 
Zolesio /]/, with the specific notation of Buscaglia and Ausas t r 10j. In particular, the surface 
gradient Vr/ of a function f : T —* R at the point x e T c R 3 is the three-dimensional 
vector 

V r / (x) = V/(x) , (5) 


where 


/( x ) = /(n r x), 


( 6 ) 
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Ilrx being the normal projection of x e M 3 onto Ik Surface gradients of vector fields are 

12 computed one Cartesian component at a time f(VrU )ij = ( Wuf)j). 

If the velocity of the surface u e V satisfies a variational formulation 

a(u,v) = /(v) VveV, (7) 

13 for some continuous bilinear and linear forms a(-,-) and /(•), then the mixed formulation 
M that enforces the inextensibility constraint is: 


Mixed inextensible formulation: Find (u,u) e VxE such that 


a(u, v) + 6(v, a) = /(v) Vv e V 

(8) 

b( u,0 = 0 VCeS, 

(9) 

where 


K v >0 = Vr • v • 

(10) 


The Lagrange multiplier a is the reaction to the inextensibility constraint, a scalar field 
that is known as surface tension. It is well known [ 3 , 11] that a necessary condition for (18]) - 
(|SD be well posed, and in particular for a to exist and be unique, is the inf-sup condition 


inf sup -— 77 ’ = a > 0 . ( 11 ) 

O^eS o^v g V ||v||x||£|| E 

We assume that the velocity space V is a subspace of 


X = {v e H 1 (T) \v n = vne H m (T)} , (12) 


15 n being the unit normal, with two possibilities: 


• If m = 1, we have X = H 1 (T) and we set 


v||x = IMIo + ^|v|i . 
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( 13 ) 
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• If m > 1 the space X has extra regularity in the normal component, and the norm is 


chosen as 


| v II x = ||v||o + ^|v|i+r*k 


(14) 


Above, || • || m (respectively, | • | m ) denotes the usual norm (respectively, seminorm) on H m {T). 
The same notation is used, without risk of confusion, for the norm and seminorm of H m (T) 
(space of vector fields with components in H’ m (Y)). The length scale i is a constant intro¬ 
duced to make the units consistent. The added regularity in the normal component results, 


in the case of lipid membrane models, from a curvature dependent energy 


12 


, i3i,y,ia 


m 


Our interest in this article lies in identifying appropriate combinations of spaces V c= X and 
£ so that is satisfied. 

Notice that b(v, a) in (jHj) is the dynamical action of a. Using the integration by parts 


formula for closed surfaces 19 


, 3 , 


16] 


a V • v = v • (— Vrcr + H a n) 
r Jr 


Vv e X, Va e H\T) , (15) 


where H = Vr ■ n is the mean curvature, one recovers the classical expression of the surface 
tension force 

f a = Vr<r — H a n . (16) 

The inextensibility constraint is frequently imposed together with the constraint of iso- 
choricity , that is, of preserving the enclosed volume. Many physical situations admit such an 
idealization, as for example the situation in which an impermeable material surface encloses 
an incompressible medium. In the case of lipid membranes, isochoricity is a consequence of 
osmotic equilibrium 8]. 

The rate of change of the enclosed volume V when T moves along the velocity field u is 

dV f - 

M = J r Un (17) 

The Lagrange multiplier that enforces this constraint turns out to be an internal uniform 
pressure p e 1 acting as a uniform normal force pn(x) at each x e T. 
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If just the isochoricity constraint is imposed, the corresponding inf-sup condition is 


inf sup 


r $ r v • n 


> 0 . 


( 18 ) 


veV |r| ||V||x 

30 If both the iuextensibility and isochoricity constraints hold simultaneously, the mixed 

3 1 formulation becomes: 


Mixed inextensible-isochoric formulation: Find (u, a, p) e V x £ x M. such that 

a ( u , v) + c(v, (cr,p)) = /(v) Vv e V (19) 

c(u, (£,g)) = 0 V (£, g) e Exl, (20) 

where 

c(v, (f,g)) = J^(£V r -v + gv-n) . ( 21 ) 


32 

33 

34 


In this case both its surface tension a (as a function of x e T) and its internal pressure 
p are uniquely defined in SxR, under the following condition on V-E, 


inf 


sup 


c(v,(f,g)) 


.I || || , . || lUf || || || . , \ | 

(?,9)sSxR veV II v||x ||(£, <?)||exR (£,g)sExR veV II v ll X || (£ , <?) IIE xR 

In what follows, we will select appropriate spaces £ and prove the stability inequalities 
(ITT]) and (1221) considering two possibilities for V: (a) the shape T is fixed, implying that 
V c X consists solely of tangential fields; and (b) the space V is unconstrained, i.e., V = X. 


$r(£ V r' v + gv-n) 


= $ > 0 . ( 22 ) 


35 2. Inf-sup conditions for purely tangential motions 

If a vector field v e X is decomposed into its tangential component v r and its normal 
part v n n, i.e., 

V = v r + v n n , (23) 

its gradient takes the form 

V r v = Vrv T + n (x) Vru„ + u n Vrn (24) 
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and its surface divergence is given by 

Vr • v = V r • v r + v n Vr • n . (25) 

The last term in (l24j) contains the curvature tensor 

H = Vrn . (26) 

36 If V only consists of tangential motions, then v n = 0 for all v 6 V, the numerator inside 

37 the inf-sup in (Il 8 p vanishes identically and thus the condition is not satisfied. The internal 

38 pressure is not uniquely defined, and modifying it has no effect on the motion of the surface 

39 continuum. 

The surface tension a, on the other hand, is well defined in L 2 (T) up to an arbitrary 
additive constant. It is uniquely defined, for example, in 

iyr). {, 6 L\T) I J), - 0 } . (27) 

Prop. 2. //V is the closed subspace of IK consisting of purely tangential motions, i.e., 


V = {v e X | v • n = 0 a.e. in T} 


(28) 


and H = Tq(T) , then the inf-sup condition Oil) holds with 


a = 


c r £ 


(29) 


40 where c r is the elliptic regularity constant (see below). 

Proof. This is proved in the same way as ([ 2 ]) is proved. Given £ arbitrary in Lq(T), let 
(p e i^ 1 (T) n Lq(T) be the solution of 

Ar<^ = £ • (30) 


Taking v = Vr^, one has 


K v >0 = ) r ^ V r ' v = Mo 


16] 


On the other hand, from the regularity estimate 

IMIo + £ Mi + ^ M2 ^ c r £ 2 Hello 

41 (where again £ is introduced to render the units consistent) one has that ||v||x 

42 £ |v|! = \(p\i + £ M 2 ^ c r £ H^llo; an d thus the claim (fTTj) is proved. 


(31) 


(32) 

Ho + 

□ 
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3. The inf-sup condition for arbitrary motions 

If a surface can move along its normal direction then condition (fIS|) is seen to hold simply 
taking v = n. The preservation of the enclosed volume uniquely defines an internal pressure. 
Turning to the inextensibility condition, one has the following: 

Prop. 3. If V = X and £ = L 2 (T), then 071) holds. 


This means that the surface tension is completely defined on an inextensible surface that 
can move freely in space. We will prove this proposition after proving Prop. U since by then 
it will become straightforward. 

Consider now the case of a surface that both is inextensible and preserves the enclosed 
volume. Rewriting the corresponding inf-sup condition (122|) as 

$r K V r • v r + (H£ + q ) v n ] 


inf 

K,9)eSx I 


sup 

veV 


Mix II(C,?)IIq> 


= p > o 


(33) 


one immediately sees that these two conditions are linearly dependent if the surface has 
constant mean curvature (a sphere). Taking q = H (constant) and £ = — 1 (constant), and 
remembering that the integral of the surface divergence of a purely tangential field over a 
closed surface vanishes (a consequence of (1T5]) taking a = 1), one has the numerator inside 
the inf-sup equal to zero for all v. For a sphere, thus, inextensibility implies the preservation 
of the enclosed volume. The internal pressure and the mean surface tension are not uniquely 
defined. The gradient of the surface tension, on the other hand, is well determined. 

For all surfaces other than spheres, however, one has: 


Prop. 4. If the surface T is not a sphere, then ( dj| ) holds with V = X and £ = L 2 (Y). 


The proof is quite straightforward, but let us for clarity first state a couple of preliminary 
facts. Let ~g denote the mean value of g e L 1 (r). 


Lemma 5. The expression 

ll(£> g)|L 2 (r)xR 


£-£« l+£ 4 (Hf + q) 2 + eti 2 


(34) 


is indeed a norm on L 2 (r) x R. □ 
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As a consequence of (l24ji one also has: 


Lemma 6. 


IMIo + ^ l v |i ^ || v t-|| o + ^ | v r1 1 + (1 + ^||H||i<»(r))||n n ||o + ^|^n|i (35) 

64 for all v e H 1 (r). 

65 Now we proceed to prove Prop. [H 

Proof. Let (£, q ) be an arbitrary element of S x R. Taking ip e iL 1 (T) n Lq(T) as the unique 
solution of A r </? = £ — £, we select v e V as 

v = Vrf + v n n , with v n = k\(H f + q) + k 2 (H — H) f . (36) 


The positive constants ki and k 2 are left unspecified for now. Rewriting 


c(v, (?,?)) = li Vr ■ V, + !//'■ - q) >;„] 


(37) 


and using (l36]i one has 


c(v, (£,?)) = U-tfo + h 


\r\(HZ + q ) 2 + (HZ + q) (ff,£-tf£) 


+ 


+k 2 


r L 


^2 


(H - H) 2 £ +(H-H)H^-0 + (H-H)(HZ + q)£ .(38) 


67 Noticing that the last term in the second integral cancels out and using that — H f) ^ 

ea ||H — H ||o ||£ — f|| 0 , one arrives at 


, t2 


c(v,(e,g)) ^ u-c\\i + h\r\(H^+ q ) 2 + k 2 \\H-H\\ 2 0 e 

-h\\H - H\\ 0 \H£ + g| ||e - ?||o ~ h\\(H - H) 2 1| 0 |£| ||£ - £||o (39) 


69 which using Young’s inequality twice yields 
c(v, (£,<?)) ^ (l-h 


H-H\\ 2 . || (iY — YY) 2 ||g 


2 1 r | 


— ko 


+ki -LJ- (H e + q) 2 + k 2 


2\\H~H\\l 
\\H~H\\ 2 ■ 


lie-as+ 


■r. 


2 


(40) 









70 


Let us now choose 


7 i gives 


h 

h 


min 


min 




c(v,(£,g)) > \u-t\\l + h^(Hti + q ) 2 + k J H 2 m h 2 
^ A ||(£,9)llL 2 (r)xR > 


(41) 

(42) 

(43) 

(44) 


72 with 


. [1 | r | 2 | r | \\ h - h \\* \\ H ~ H \\ 2 0 ] 

mm { 2 ’ 4 t 4 \\H ~ H \\ 2 0 ' 2 P ’ 4 £ 2 \\(H ^ H) 2 \\f 2 j 

73 At the same time, from (j35|) and the estimates 

||v T ||o + ^|v T |l < C r £||£-el0 


(45) 


(46) 


u n || 0 < e\V\^\Hi + q\ + P\\H-H\\o\Z\ (47) 


K\k < t 2 \H\ k \Z\ 1 (48) 

74 it follows that 

i|v||x < Cr£\\€ - €||o + (1 + £ ||H|| L co (r) ) t |r|5 \HZ + q\ + 

+ [(l+£ ||H|| L oo (r) ) e 2 1 \H - H\\ 0 + e \H\, + r +2 1 H\ m ] lei 
^ B ||(^?)||L2 ( r)xR (49) 


with 

B 2 = cj e + (1 + l || H | n « (r) ) 2 | r |+ [(1 + £ IIHIU.jr,) l \\H - ff || 0 + P\H |, + r +1 |ff |„] 2 

(50) 

75 The claim is proved with /3 = A/B. □ 
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Notice that A is equal to zero for a sphere, and thus j3 = 0. However, Prop. [3] is true 
irrespective of Y being a sphere or not. Let us modify the previous proof to prove it. 


Proof, (of Prop. [3]) Taking the same v as before, and particularizing 


to q = 0 one gets 


&( v >0 > 1 kl 




-2 


? > cillli 


(51) 


with 


C = min 


1 ki H 2 k 2 \\H-H\\ 2 0 

2 ’ 2 1 2 in 


and now C > 0 even if H = H. Proposition [3] is thus proved with £ 


L 2 (r) and a 


(52) 

C/B. 

□ 


4. Discrete inf-sup condition 


The continuous inf-sup conditions proved above also allow for the extension to surface 
finite elements of the discrete counterpart known as Verfiirth’s lemma [l7]. It is central 


in the numerical ana’ 


Squares methoc 


method 

[ 5 , 3 . 






ysis of stabilized finite element methods such as the Galerkin-Least- 


19 ], A variant of one such method, the pressure gradient projection 


22 j, has recently been successfully implemented for lipid membrane models 


Let V/j c V and c £ be surface finite element spaces, as defined in Dziuk and Elliott 
[u]|. Notice that these are lifted spaces, which are defined with the aid of a faceted surface 
but consist of scalar functions (or vector fields) defined on the “exact” surface T. 


Prop. 7. If the space T,h consists of continuous functions, then there exist 7 > 0 and S > 0 , 
independent of the mesh size h = max^ hx, such that 


sup 

v h eV h 


$ r (^ V r • Vfe + gv/i • n) 
II v ft. I x 


^ 7 ||(6i, <?)|sxR 



v r^||£2 (K) 


(53) 


for all (£ ft , q) e £ h x R. 
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The proof assumes the existence of a Clement-type quasi-interpolatory operator Zh : 


V —*• V/ t satisfying 


[23, \n\. 


Mx < Cl |jv||x (54) 

||v -X h v\\ L2{K) ^ c 2 Hk || V r v||i 2 (a;jr) (55) 

where uJ k is the union of all elements that share at least one node with element K 

Proof. For arbitrary Ci and q, from Prop. Owe have 

awti- Ml / c(Z h v, (£ h ,q)) , c(v - Z h v, (£ h , q)) 

p ||Ufe,g)||ExR < sup-jj—r;- + sup-jj—jj- ■ (56) 

veV || v || X veV || v ||x 

It is clear that 

c(X h v : (f h ,q)) . c(l h v,(£ h ,q)) c(v h ,(£h,q)) 

sup--— -^ ci sup-—— -^ ci sup -- —jj- . (57) 

veV ||v||x vgV ||2^V|| X v k eV h ll v h||x 

Concerning the second term in the right-hand side of (l56]h and denoting w = v — I^v, we 

have 


c(v-Z h v,(£ h ,q)) = 


Y [f h Vr-w + qHw-h] 

'cT JK 


K sTh 


Y Ci w • iy T V (-VrC • w + fhHw ■ n) 

' C T 7 JdK KeTh 


(58) 


K <- T,, 


The first sum vanishes if £/i is continuous across element boundaries and T is C 1 . As a 
consequence, since ||wj| L 2 (x) ^ c 2 h K || V r v||L 2 (a> K ), 


c(v-X h v,{£ h ,q)) ^ j 


Y ^|Vr^||L W ) + h \\H ||n°°(r) Iloilo 

\K e Th 


Mix (59) 


where 


C 3 = c 2 1 sup 

vgV 


TjKeT h 

II v l|x 


Replacing into (j55j) one gets 


own mi ^ c(v h ,(£ h ,q)) c 3 

p||(^h,g)||sxR < ci sup -j :—jj- + — 

v h eV h ||Vh||x t 


D PjVr&ll 

\K e Th 


L 2 (K) 


(60) 


+ ZiEZnn | & | 0 


£ 


( 61 ) 


which proves (153]) taking 7 = j3/{2 ci), 5 = c 3 /(£ci) and h small enough. 


□ 
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5. Concluding remarks 


It has been shown that the inextensibility of a surface continuum, analogous to the 
incompressibility of a volumetric medium, is a well-posed constraint for closed surfaces 
evolving in M 3 . It gives rise to a surface tension field a that is uniquely defined in L 2 (T). 

The simultaneous imposition of both the inextensibility constraint and the isochoricity 
constraint (preservation of the enclosed volume) is also well posed, the only exception being 
that of a spherical configuration of the surface continuum. In all other cases, both the 
surface tension field a e L 2 (T) and the internal pressure p e 1 are uniquely determined by 
the constrained problem. 

The estimates in the proofs require the surface T to be of class C 2 and have its mean 
curvature H in H m (T), m ^ 1. 

On the basis of the exact well-posedness, a discrete stability result was established for 
discretizations of the surface tension field that consist of continuous interpolants. It consists 
of a modified inf-sup condition (sometimes called VerfiirtlTs lemma) which plays a central 
role in the justification of stabilized methods for incompressible flow. The extension of this 
inf-sup condition to deformable surfaces thus justifies the stabilized treatment of the surface 
tension proposed in recent work on lipid membranes 0 , 01 - 
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